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ABSTRACT 

The nebular spectra of the energetic Type Ic supernova SN 1998bw 
(hypernova) are studied. The transition to the nebular phase occurred at an 
epoch of about 100 days after outburst, which is assumed to coincide with 
GRB980425. Early in the nebular epoch the spectra show the characteristics of 
a typical SN Ic spectrum, with strong lines of [O i], Ca II and Mg i], and lines of 
[Fe n]. However, the [Fe n] lines are unusually strong for a SN Ic. Also, lines of 
different elements have different widths, indicating different expansion velocities. 
In particular, iron appears to expand more rapidly than oxygen. Furthermore, 
the [O i] nebular lines decline more slowly than the [Fe n] ones, signalling 
deposition of 7-rays in a slowly-moving O-dominated region. These facts suggest 
that the explosion was aspherical. The absence of [Fe ill] nebular lines can be 
understood if the ejecta are significantly clumped. A schematic picture of what 
this very unusual stellar explosion may have looked like is presented. 

Subject headings: supernovae: individual: SN 1998bw - nucleosynthesis 
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1. Introduction 

SN 1998bw was a very interesting object from the word go (e.g. Nomoto et al. 2001). 
Discovered in the error box of GRB980425, and very possibly linked to it, this object 
was soon identified to be a supernova (SN) from its light curve, which was different from 
those of typical optical transients of GRB's. Extensive data were collected at ESO, where 
SN 1998bw was designated a Target of Opportunity (Galama et al. 1998; Patat et al. 2001), 
and at other observatories. 

Early spectra were rather blue and featureless, quite unlike those of other known SNe. 
However, a somewhat more careful inspection showed clear similarities with the spectra of 
Type Ic SNe (SNe Ic), but with one major difference: the absorption lines were so broad in 
SN 1998bw that they blended together, giving rise to what could even be confused with an 
emission spectrum. Very large expansion velocities (~ 30, 000 km s _1 ) were measured in 
the ejecta. Also, SN 1998bw was very bright for a SN Ic: for a redshift distance of 39 Mpc 
(z = 0.0085 and H = 65 km s _1 Mpc -1 ), the SN reached a maximum V= —19.4 mag, 
which is bright even for a SN la, while a typical SN Ic like SN 1994 1 was almost 2 mag 
dimmer (e.g. Nomoto et al. 1994; Richmond et al. 1996). 

Modelling of the early data (Iwamoto et al. 1998, hereafter IMN98) led to several 
striking conclusions. First, the SN produced about 0.7M & of 56 Ni, as much as a SN la. This 
was necessary to power the light curve via the deposition of the 7-rays emitted in the decay 
of 56 Ni into 56 Co and 56 Fe. 

IMN98 discuss that the light curves are degenerate, i.e. different combinations of ejecta 
mass M e j and kinetic energy E kin can all reproduce the data equally well. Spectra came to 
the rescue then: spectrum synthesis showed that the only way to get the observed extensive 
line blending was to have significant amounts of material moving at very large velocities. 
Since models with different E kiD give different degrees of blending, it was possible to select 
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as best model one with M ej = 1O.9M and E kin = 3 • 10 52 erg. This had an exploding core 
mass of 13.8M and was designated C0138 accordingly. A typical SN Ic has M e - S = 1M , 
and the typical kinetic energy of most SNe of all types is only about 10 51 erg (aka 1 foe). An 
even higher value of 50 foe was obtained for SN 1998bw by Nakamura et al. (2001a). This 
gives a better fit to the declining part of the light curve and to the spectra. A comparable 
value was obtained by Branch (2001) from an analysis of the early time spectra. Because of 
its exceptionally large -Ekin, SN 1998bw was called a 'hypernova'. The progenitor must have 
been a massive star, and we estimated a main sequence mass of about 4OM . Also, the 
remnant mass, which is computed by allowing only as much of the 56 Ni produced by core Si 
burning to be ejected as is necessary to power the SN light curve, turned out to be ~ 3M 
Nakamura et al. (2001a,b). This exceeds the maximum mass of a Neutron Star, suggesting 
that the explosion that was observed as SN 1998bw resulted in the creation of a Black Hole. 

Since SN 1998bw was probably connected to a highly non-spherical event like a GRB, 
departure from spherical symmetry could be expected. Early polarization measurements 
confirmed this. Polarization of ~ 0.5%, possibly decreasing with time, was detected (Kay 
et al. 1998, IMN98, Patat et al. 2001). This was interpreted as an axis ratio of about 2:1 in 
the expanding ejecta (Hoflich et al. 1999). Models of the collapse of a rotating stripped core 
of a massive star were developed, which confirmed that the explosion should be asymmetric, 
thus establishing a possible link with the GRB (MacFadyen & Woosley 1999). 

All the above conclusions were based on only about two months worth of data. Patat 
et al. (2001) presented extensive data covering more than one year. Note that here and 
throughout the paper epochs refer to the SN time since outburst, which is taken to have 
coincided with GRB980425, corrected for time dilation (the host galaxy of SN 1998bw, 
ESO 184-G82, has a redshift cz = 2532 km s" 1 ), while Patat et al. (2001) used the time 
from B maximum, which occurred on 10 May, 1998). The evolution of SN 1998bw after 
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the first two months reserved several surprises: after day 60 the light curve began declining 
less steeply than the model prediction (Mc Kenzie & Schaefer 1999); the nebular spectrum, 
which started to develop after about day 100, initially appeared to be a composite of a 
SN la and a SN Ib/c spectrum, showing strong Fe lines in the blue, typical of SNe la, and 
strong O and Ca lines in the red, typical of SNe Ib/c. What was even more peculiar was 
that there were two families of lines: one with a large velocity comprising mostly the Fe II 
lines, and another with a smaller velocity ([O i], Mg rj). Also, the various lines declined 
with different rates, which apparently contradicts the hypothesis that the only heating 
source is the decay of 56 Co (Patat et al. 2001). 

Several papers have addressed the behaviour of the light curve over the first 200 
days (Mc Kenzie & Schaefer 1999), 600 days (Sollerman et al. 2000) and up to 800 days 
(Nakamura et al. 2001a). In all these papers the suggestion has been made that either the 
density distribution computed in spherically symmetric explosion models is not correct or 
that some degree of asymmetry in the explosion is probably required to explain the unusual 
behaviour of the light curve. In this paper we focus on the late time spectra, and show that 
they too require significant asymmetry in the ejecta. 



2. Nebular spectra 

As discussed by Patat et al. (2001), the transition from photospheric to nebular spectra 
took place between about day 65 and day 115, although net emission may have been present 
in the Ca II IR triplet as early as day 37. Later other Ca II lines, the H&K doublet and the 
forbidden 7320A line became visible, followed by lines of [O i], Mg i], by Na I D and by an 
emission at ~ 5200A. This is most likely due to a complex of [Fe n] lines, although at the 
earliest phases a possible contribution is present from permitted Fe II transitions, which 
can produce nebular emission at sufficiently high density (Filippenko 1989), and most likely 
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also from [O i] 5577A. At about day 115 the spectrum of SN 1998bw was an interesting 
hybrid between that of a SN la (strong [Fe II]) and that of a SN Ib/c (very strong [O i], 
Mg i], Ca n). This is already an indication that both Fe and the a-elements are present 
with significant abundances in the ejecta. Another SN Ic for which an intense blue flux has 
been observed in the nebular epoch is SN 1999dq (Matheson et al. 2000a). This blue flux 
has been attributed to Fe II lines, but the wavelength coincidence of this feature with the 
[Fe n] one in SN 1998bw makes the [Fe n] identification also a serious possibility. SN 1999cq 
was a bright SN Ic (at least using a redshift distance), so the presence of a larger amount of 
Fe than in 'normal' SNe Ic is to be expected. 

The nebular spectra span epochs from 108 to 388 days. In all the spectra six dominant 
features can be identified: Mg i] at 4500 A (1), the [Fe il] complex near 5250A (2), Na I D 
at 5900A (3), [O i] 6300A (4), a [Ca n]-[0 il] blend near 7200A (5) and the Ca II IR triplet 
near 8500A (6). However, although these features persist for almost one year, the change in 
the properties of the spectrum are striking. In particular, the [O i] and Mg i] lines become 
stronger relative to the [Fe il] lines as time goes on, but at the same time they also become 
narrower. While in the first spectrum all lines are broad (v ~ 12000 km s _1 ), already in the 
second spectrum [O i] , Mg i] and partly also Ca II IR are narrower than the [Fe il] complex. 
This behaviour continues for the duration of our spectral coverage, throughout which the 
[O i] and Mg i] lines and the [Ca n]-[0 il] blend grow in strength relative to the other 
three features. In practice, the decline rates of the various emission lines are different, so 
that while in the first spectrum features (1), (2), (4) and (6) all have similar peak fluxes, 
and feature (5) is as strong as feature (3), in later spectra [O i] is the strongest line by 
far, while the [Fe il] feature rapidly decreases in strength relative to the other lines, so 
much so that after about one year the spectrum of SN 1998bw would be very difficult to 
distinguish from that of a normal SN Ib/c (see Patat et al. 2001 for a direct comparison and 
Matheson et al. 2001 for lots of late-time spectra of SNe Ic). At late epochs the strength 
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of the [O i] and Mg i] lines is mostly due to the narrow component, which develops with 
time but is absent in the [Fe il] feature. All the lines that decline slowly become eventually 
significantly narrower than the [Fe n] features. The weak Na I D line also has a narrow 
profile, and declines relatively slowly, while the Ca II IR triplet is almost as broad as the 
[Fe il] feature and declines almost as fast. The decline of the ratio Ca II IR/[Ca n] 7300A is 
to be expected in a nebula of progressively decreasing density (Ferland & Persson 1989). 

The features described above persist for the entire period sampled, the only noticeable 
change being that the [Fe il] feature moves bluewards by about 130A between day 108 and 
day 214. This may be due to a decreasing contribution of an underlying 'photospheric' 
component, which is probably still significant at blue wavelengths in the earlier epochs, 
and to a progressively decreasing contribution of permitted Fe II emission, as the density 
decreases, and of [O i] 5577A, which has a rather high excitation potential (its lower level 
is the uppper level of [O i] 6300A), and whose strength should decline as the electron 
temperature drops. After day 214, though, the emission feature stops moving, and its 
identification as the same [Fe il] feature observed in SNe la is beyond doubt (Patat et al. 
2001, Fig. 11). However, the decline rate of the luminosity of the [Fe il] feature remains 
larger than that of both [O i] and Mg i] at all epochs (Patat et al. 2001, Fig. 14), and so 
this is not just due to changing contributions to the [Fe il] feature. 

Another feature which receives different contributions is Na I D. Early in the nebular 
phase the ground state [Co ill] multiplet near 5900A contributes to this line, as it does in 
SNe la. However, the contribution is smaller in SN 1998bw because of the lower degree 
of ionisation. In fact, in our models Na I D is responsible for about 80% of the observed 
feature. At later times the contribution of Co ill is further reduced because 56 Co decays to 
56 Fe. This feature may also receive a contribution from He I 5876A, which may be excited 
non-thermally since He and 56 Ni coexist in the ejecta (Maeda et al. 2001, Fig.l). 
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At all epochs, the strength of the [Fe n] feature is unusually large for a SN Ib/c. 
Patat et al. 2001 (Fig. 11) compared the nebular spectrum of SN 1998bw with that of the 
SN Ic 1996aq at an epoch of about 9 months, showing that in SN 1998bw the entire blue 
part of the spectrum, where the [Fe n] lines are important, is much stronger relative to 
the [O i] and Ca II lines in the red. Spectra of SNe Ic at an epoch of about 5 months, 
corresponding to the epoch when the spectrum of SN 1998bw had just turned nebular, have 
been published by A. Filippenko and collaborators (Filippenko et al. 1990 (SNe 1985F, 
1987M); 1995 (SN 19941); Matheson et al. 2001 (SNe 1990U, 1990aa)), and a comparison 
by eye of their spectra to those of SN 1998bw confirms the relative strength of the [Fe II] 
feature in SN 1998bw even at the earlier epoch. 

These properties of the late-time spectrum of SN 1998bw are odd in many respects. 
First, all (spherically symmetric) evolution-collapse-explosion models predict that O, which 
is found in the unburned part of the ejecta, is located above Fe in mass coordinates. 
Therefore O should have a larger velocity than Fe if anything, but not a smaller one. 
Secondly, if the spectrum is powered by 7-ray deposition, the region of the ejecta which 
should feel the heating the most should coexist with Fe, which is mostly the product of 
the decay of 56 Co. Therefore, all lines should have a similar width, that of the Fe lines, 
and they should decline at about the same rate. This is the case for SNe la. The next odd 
fact is that SN 1998bw had broad [Fe 11] lines like SNe la, but lacked the strong [Fe ill] 
peak near 4700A, which is the strongest line in the nebular spectra of SNe la. If the [Fe 11] 
identification is correct, the absence of [Fe ill] clearly must be the result of a different 
(lower) degree of ionisation in SN 1998bw. The above facts indicate that something unusual 
must be going on with the way the spectrum is powered. 

We used a NLTE code to compute the nebular spectrum of SN 1998bw. The code 
computes the deposition of the 7-rays from the decay of 56 Co to 56 Fe in a nebula of uniform 
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density and composition. Together with the positrons, which are also produced in the 
decay of 56 Co and which are assumed to be trapped, this provides heating to the nebula. 
Cooling by nebular line emission is computed, a matrix of level populations is solved in 
NLTE and an emission spectrum is obtained. This approximation is reasonable for SN la 
nebulae (Axelrod 1980), so we apply it also to SN 1998bw. It is reasonable to do so as long 
as the nebular spectrum of SN 1998bw is mostly powered by 56 Co decay. Furthermore, the 
density in the explosion model C0138, which reproduces the early light curve and spectra 
of SN 1998bw, is actually flat below about 12,000 km s" 1 (Nakamura et al. 2001a), and 
mixing in the ejecta appears to be required because a stratified abundance distribution 
gives nebular lines with a flat-top profile (see Mihalas 1978, p. 477), which is inconsistent 
with the observations (Sollerman et al. 2000). Mixing was also adopted in the light curve 
models of IMN98. Masses of the various elements are given as input to the code. 

We modelled 5 spectra of SN 1998bw, dating from three months after maximum, when 
the SN was still entering the nebular phase, to more than one year after maximum. 

Since our code assumes a uniform density, it can be used to determine the outer 
velocity of the nebula by fitting the width of the emission features, including complex 
blends like [Fe il], which is the broadest feature in all the nebular spectra of SN 1998bw. 
Because line blending is computed, the comparison of the widths of synthetic and observed 
emission features is realistic. We found that the velocity of the [Fe il] feature decreased 
only marginally, from 12,000 km s" 1 in the first spectrum to 10,000 km s" 1 in the last two, 
indicating that 7-ray deposition slowly became less efficient in the outer part of the nebula. 

The major difficulties in reproducing the spectra were the presence of the narrow lines, 
in particular [O 1] and Mg 1], and the apparent absence of [Fe ill] lines. As for the narrow 
lines, we chose to fit their peak flux, thus deriving an upper limit to the masses of these 
elements. The absence of [Fe ill] lines is interesting. Typical models of the nebula have 
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electron density n e ~ 10 7 cm -3 , electron temperature T e ~ 7000 K, and a ratio Fe Ill/Fe II 
~ 0.5, so that both [Fe n] and [Fe ill] lines are strong. This situation is typical also in 
SNe la. The only way to decrease the strength of the [Fe III] lines without also reducing 
that of the [Fe n] lines is of course to decrease the ionisation. This can only be achieved 
by increasing n e to favour recombination. Iron is not the main contributor to the electron 
density because of its relatively low abundance, so increasing n e does not mean increasing 
Fe Ill/Fe II significantly. n e cannot be increased by simply increasing the density, as this 
would imply increasing the mass, which would result in more emission in all lines. 

One solution is to introduce clumping in the ejecta. This is done using the concept of 
volume filling factor. In practice, the state for a mass M is computed using a volume fV, 
where V is the entire volume of the nebula, computed from the outer velocity and the time 
since outburst, and / is the volume filling factor, with < / < 1. With this approximation, 
the density of the gas is increased by a factor / to simulate clumping. The full volume V 
is then used to distribute line emissivity and to compute the spectrum, which corresponds 
to the assumption of uniformly distributed clumps. We found that for a filling factor of 
about 0.1, n e (averaged over the whole nebula) increases to ~ 3 • 10 7 cm~ 3 , while T e drops 
to ~ 6000 K, producing a lower Fe Il/Fe ill ratio of 0.05-0.1, which successfully eliminates 
the [Fe ill] lines, especially at the more advanced epochs. The presence of clumping is 
not surprising in the ejecta of core-collapse SNe (e.g. Spyromilio et al. 1993 for SN 1987A 
and Spyromilio 1994, Matheson et al. 2000b for SN 1993 J). Sollerman et al. (2000) using 
a detailed density distribution in spherical symmetry obtained a lower degree of ionisation 
in the Fe-rich part of the ejecta, but they had to reduce the ejecta velocity artificially with 
respect to the results of the explosion models in order to obtain a good agreement with the 
width of the observed lines. Our models for the epoch they also modelled, day 139, are 
similar to theirs, providing perhaps a slightly better fit in the blue region where the [Fe ill] 
lines form. We also confirm the absence of lines of Fe I. 
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The properties of our models fitting the broad [Fe n] feature are summarised in Table 
1, and the synthetic spectra are compared to the observed ones in Figure 1. 

For all the 'broad' models, the 56 Ni mass was about O.65M . This is in good agreement 
with the light curve calculations, although it is perhaps a slight overestimate since there 
is more flux in the synthetic spectra than there is in the observed ones, owing to the fact 
that the narrow observed lines are reproduced as broad lines. It is also encouraging that 
the 56 Ni mass is very nearly constant over the entire period, showing that the [Fe n] feature 
is indeed due to 7-ray deposition in an expanding nebula. 

The O mass, on the other hand, ranges between 3 and 5M in the various models, 
increasing first and then decreasing. In the first spectrum the profile of [O i] 6300A is 
broad, and so our estimate (M(O) ~ 3M ) is probably more realistic, considering the 
large errors which are always incurred in when determining the mass of O from that line 
(Schlegel & Kirshner 1989). However, at later times the narrow component dominates, our 
broad synthetic [O 1] lines have a larger EW than the observed lines, and our O masses are 
certainly overestimated. At very late phases the O mass drops again to about 3M Q , but 
this may be just an artifact of using a constant density model with a given velocity radius 
when the real ejecta had a decreasing density and so effectively a velocity radius decreasing 
with time. 

The fact that the narrow lines of [O 1] and Mg 1] decline less rapidly in luminosity 
than the [Fe 11] feature with time appears to indicate that there is a significant amount of 
these elements at low velocity. The ratio of the [Fe 11] feature and the broad component 
of [O 1] 6300A in the spectrum at day 108 also suggests that the Fe/O ratio is larger at 
high velocity. All this is very difficult to reconcile with the classical picture of a spherically 
symmetric explosion, even taking into account mixing. 

Sollerman et al. (2000) computed a nebular spectrum for day 139 using the original 



- 12 - 



density and abundance distributions of model C0138 (IMN98), and obtained broad [01] 
lines with a flat emission top. This is a consequence of distributing O in a shell at high 
velocity. They overcame this problem by mixing the composition to eliminate the flat tops 
and by arbitrarily decreasing the velocity of the outer part of the nebula to reduce the 
width of the [O i] lines. However, our code can only handle nebulae of uniform density. 
Thus, in an alternative set of models we tried to fit the narrow-line spectrum, assuming 
it is also powered by 7-ray deposition. Such models were computed only for the last four 
epochs, since in the first spectrum the narrow [O 1] line is still weak. The parameters of the 
models are detailed in Table 2, and the synthetic spectra are shown in Figure 2. 

The synthetic narrow-line spectra fit well the lines of Mg 1], Na I D and [O 1]. The 
width of the lines decreases only slowly with time, as was the case for the broad [Fe 11] 
feature: velocities range from 7500 to 5000 km s -1 . The synthetic [Fe 11] feature is in these 
models significantly too narrow. Therefore we chose to reproduce its peak only. Once 
again we found that the 56 Ni mass is very nearly constant in all four models, with a value 
of ~ 0.35M Q . However, this is certainly an underestimate, as the total flux in the [Fe 11] 
feature is not reproduced. Values of the O mass in the low velocity nebula range between 
1.5 and 2.7M , with the estimate decreasing with time after day +201. O masses which 
decrease with time were also derived by Schlegel & Kirshner (1989) for the SNe lb 1984L 
and 1985F and for the SN II 1980K. 

An upper limit to the O mass can also be derived from the flux in the [O 1] 6300A line 
and the nebular temperature in the high density limit (n e > 10 6 cm -3 ) using the formula 

M(O) = 10 8 d 2 F([O l])e 2 - 28/T4 (1) 

(Uomoto 1986), where M(O) is the O mass in M , d is the distance to the SN in Mpc, F is 
the flux of the [O 1] line in erg cm~ 2 s _1 and T 4 is the temperature of the nebula in units 
of 10 4 K. The high density limit should apply given the conditions in our model nebulae 
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(Tables 1 and 2). 

The values we obtained for SN 1998bw (d = 39.3Mpc) using values of T 4 as derived 
from our models are listed in Table 3. The O mass thus derived is in good agreement with 
that used in the model calculations, except at the latest epochs, when the density is lower 
and the condition of Eq.(l) is not as well satisfied. 

3. A 'portrait' of SN 1998bw 

Most of the work on SN 1998bw has been based on spherically symmetric models. The 
very bright and relatively broad light curve of SN 1998bw can be reproduced by a family 
of explosion models with various values of the fundamental parameters (M e j, £idn)- All 
models require M( 56 Ni)~ 0.4 — O.7M to power the bright light curve peak. This is about 
one order of magnitude larger than in typical core-collapse SNe. Models with different E kin 
yield different synthetic spectra, and by comparing with the observed early-time spectra of 
SN 1998bw and trying to fit the very broad absorption features, Nakamura et al. (2001a) 
selected a model with M cj = 1O.9M , E kin = 5 • 10 52 erg (model CO138E50). This model has 
the same M cj but a larger E kiQ than the model CO 138 of IMN98, and it yields a better fit 
to the observed bolometric light curve, the evolution of the photospheric velocity and the 
spectra of SN 1998bw. The large value of E kin easily qualifies SN 1998bw as 'the' Type Ic 
Hypernova. The mass of the exploding CO star was 13.8M , which implies a main sequence 
mass of ~ 4OM . Similar results were also obtained by Branch (2001) from models of the 
early time spectra. 

These explosion models work well around maximum, but the predicted light curve 
tail is quite a bit steeper than the observed one, starting at about 60 days and continuing 
until at least day 200. This signals that either these energetic models underestimate 7-ray 
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deposition at those phases or that an additional source of energy is present, or both. 

The nebular spectra of SN 1998bw also had a peculiar evolution. Early on, the SN 
showed a 'composite' spectrum: [Fe n] lines, typical of SNe la, were strong, and so were 
lines of [O i] and Mg i], which are typical of SNe Ib/c. At the same time [Fe ill] lines, also 
typical of SNe la, were absent. The [O i] and Mg i] lines grew stronger with time relative to 
the [Fe n] lines, but this was due mostly to a narrow component, which became more and 
more the dominant one as time progressed. The emergence of the narrow profiles occured 
at about the same time as the light curve deviated from the model prediction. 

The broad spectrum could be reproduced with M ( 56 Ni)~ 0.65M Q and an outer nebula 
velocity of ~ 12,000 km s" 1 . This is consistent with the distribution of 56 Ni in CO138E50 
of Nakamura et al. (2001a), but the estimate of the 56 Ni mass is somewhat larger here. 
Note, however, that the light curve obtained with CO138E50 underestimates the late time 
flux (Nakamura et al. 2001a, Fig.5), and so this result is not unexpected. The narrow lines 
have a width of only ~ 6, 000 km s -1 , and can be reproduced with a 56 Ni mass of O.35M . 
The velocity of the narrow [O i] line appears to be consistent with the lower energy model 
CO138E10 of Nakamura et al. (2001a). 

On average, our nebular models required enclosed masses of ~ 3M below a velocity 
of 6000 km s" 1 and of ~ 5M below 12000 km s -1 . Both these values are larger than the 
corresponding values from model C0138 by about 2M Q . This may suggest the presence of 
a high density inner part of the ejecta, which could explain the slow decline of the light 
curve between days 60 and 200, as suggested by Nakamura et al. (2001a) and Chugai 
(2000). A similar suggestion was made by Iwamoto et al. (1999) and Mazzali et al. (2000) 
for the other Type Ic hypernova, SN 1997ef. However, the two situations are different. For 
SN 1997ef, evidence of a high density region at low velocities came from the fact that the 
spectra remain photospheric for a long time, and that the photospheric velocity drops below 
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the velocity of the inner cutoff. In a 1-dimensional explosion model, this is the separation 
between the material which is ejected by the SN and the infalling matter which forms the 
compact remnant. In SN 1998bw the photospheric phase ends earlier, and the need for the 
central high density region comes from the behaviour of the light curve, and it is therefore 
a less direct result. 

Rather strong clumping had to be introduced to increase recombination and keep 
the [Fe ill] lines from forming. Introducing clumping may slow the light curve somewhat. 
The increasing strength of the low-velocity intermediate-mass element emission, however, 
indicates that this will not be sufficient to explain all observations. The same would apply 
to the approach adopted by Sollerman et al. (2000), who reduced the velocity of the ejecta, 
thus enhancing 7-ray deposition, and could thus improve the fit to both the light curve and 
the spectrum. Also, with this approach the [Fe 11] lines are not expected ever to become 
broader than the [O 1] line. 

A more complete solution which takes into account all the observed facts should 
probably include the presence of asymmetry in the explosion, both geometrical and in the 
density distribution. If SN 1998bw was actually linked to GRB980425, the fact that the 
explosion was aymmetric should not be surprising. Observations indicate that a rather 
large mass of O-dominated material is present at low velocity. Spherical explosions do not 
allow that: they are very effective at 'emptying' the central region, and always place the 
unburned elements at the top of the ejecta. 

A 2-dimensional explosion model was computed by Maeda et al. (2001). If the explosion 
is highly asymmetric, it produces large quantities of unburned material (mostly C and O) 
expanding at low velocity in directions away from the axis along which most of the energy 
was released and 56 Ni synthesised. Our vantage point must have been very close to that 
axis, because we also detected the GRB, as Maeda et al. (2001) show from an analysis of 
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the line profiles. At early times, the fast-expanding lobes must have been much brighter 
than the rest, and so we observed the broad-lined spectra and the bright light curve. The 
fast-moving regions rapidly became thin, and soon emission lines appeared. Initially the 
emission lines were broad, dominated by the hyper-energetic lobes. The distribution of 
density and abundances in a 2-dimensional model (Maeda et al. 2001, Fig. 2) shows that in 
the region of the ejecta with velocities of about 10000 km s^ 1 the density contours deviate 
less from spherical symmetry. Therefore, early in the nebular phase, when these fast-moving 
regions contribute to the nebular spectrum, a sufficient amount of high velocity oxygen is 
present that the [O i] line can be expected to be broad. Later, though, the 7-rays from 
the fast-moving 56 Co could escape that region more and more easily, and the inner, slowly 
moving, O-dominated part of the ejecta would dominate the spectrum. A smaller amount 
of 56 Co is present at low velocities, but the 7-rays it produces are trapped more efficiently 
and excite the O and Mg there. At lower velocities, though, the distribution of oxygen 
becomes preferentially equatorial, and one expects the [O 1] line to decrease significantly in 
width as the effective size of the emitting nebula decreases. Maeda et al. (2001) suggest 
that the relative width of the [Fe 11] complex and the [O 1] line can be explained by such a 
model. However, because O dominates over Fe at low velocity, the strength of the [O 1] line 
should increase relative to that of the [Fe 11] feature, as is indeed observed. 

The low- velocity cutoff of the asymmetric explosion model (Maeda et al. 2001, Fig. 2) 
is only about 1500 km s _1 , which is significantly less than in the ID models (5000 km s _1 
in CO138E50; Nakamura et al. 2001a, Fig. 7), also because of the lower overall E kin of the 
asymmetric models (see below). This feature of the 2D models seems to match nicely the 
need for a low-velocity, high-density region resulting from the behaviour of the light curve. 
Additional deposition in the low-velocity region of the 7-rays produced in the high velocity 
region may in turn increase the deposition function above what our spherically symmetric 
models could estimate, and may help explain the slowly declining tail of the light curve. 
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That the SN 1998bw explosion was asymmetric is not a new suggestion: first the 
polarization measurements indicated an axial ratio of 2-3:1 (Hoflich et al. 1999), and the 
calculation of the explosion of a rotating core (McFadyen & Woosley 1999) also gave similar 
results, in an effort to explain the connection between SN 1998bw and GRB980425. More 
detailed results have to await detailed numerical models in two dimensions. If the explosion 
was asymmetric, our results for E^ in in spherical symmetry are overestimated, because those 
would only refer to the fast-moving part of the ejecta. In fact, Maeda et al. (2001) propose 
a model with E kin = 10 foe. However, the estimate of the 56 Ni mass in the nebular epoch 
(~ O.5M ) is probably correct, as it should not be significantly influenced by the asymmetry 
as long as the ejecta are completely optically thin. In this case, in fact, every optical photon 
created by the deposition of 7-rays and positrons escapes the nebula immediately and can 
be observed as SN light. As for the value of M ej , this can only be determined by computing 
light curves and spectra using 3D hydrodynamical models of the explosion (e.g. Khokhlov 
et al. 1999), but a careful analysis of the spectra, especially at late times, when both the 
fast and the slow components are observable, can yield at least approximate results. 

Acknowledgements. This work has been supported in part by the Grant-in-Aid for 
Scientific Research (12640233) and COE research (07CE2002) of the Japanese Ministry of 
Education, Science, Culture, and Sports in Japan. It is a pleasure to thank Alex Filippenko 
for enlightening discussions on the late-time spectra of SNe Ic. 



-18- 

REFERENCES 

Axelrod, T., 1980, Ph.D. Thesis, Univ. of California, Santa Cruz 

Branch, D., 2001, in Supernovae and Gamma Ray Bursts, ed. M. Livio et al. (Cambridge: 
Cambridge University Press), in press 

Chugai, N.N., 2000, Astron. Letters 26, 797 

Ferland, G.J., & Persson, S.E., 1989, ApJ 347, 656 

Filippenko, A.V., 1989, AJ 97, 726 

Filippenko, A.V., Porter, A.C., & Sargent, W.L.W., 1990, AJ 100, 1575 

Filippenko, A.V., et al, 1995, ApJ 450, Lll 

Galama, T. J., et al, 1998, Nature 395, 670 

Hoflich, P., Wheeler, J.C., & Wang, L.-F., 1999, ApJ 521, 179 

Iwamoto, K., et al., 1998, Nature 395, 672 

Iwamoto, K., et al., 1999, ApJ 534, 660 

Kay, L.E., Halpern, J. P., Leighly, K.M., Heathcote, S., Magalhaes, A.M., & Filippenko, A.V., 
1998, IAU Circ. 6969 

Khokhlov, A.M., Hoflich, P., Oran, E.S., Wheeler, J.C., Wang, L.-F., Chtchelkanova, & 
A.Yu., 1999, ApJ 524, L107 

MacFadyen, A.I. & Woosley, S.E. 1999, ApJ 524, 262 

McKenzie, E.H., & Schaefer, B.E. 1999, PASP 111, 964 



-19- 

Matheson, T., Filippenko, A.V., Chornock, R., Leonard, D.C., & Li, W., 2000a, AJ 119, 
2303 

Matheson, T., Filippenko, A.V., Ho, L.C., Barth, A.J., & Leonard, D.C., 2000b, AJ 120, 
1499 

Matheson, T., Filippenko, A.V., Leonard, D.C., & Shields, J.C., 2001, AJ 121, 1648 

Maeda, K., Nakamura, T., Nomoto, K., Mazzali, P. A., & Hachisu, I., 2001, ApJ, submitted 
( Fstro-ph/001 1003|) 

Mazzali, P. A., Iwamoto, K., & Nomoto, K., 2000, ApJ 545, 407 

Mihalas, D., 1979, Stellar Atmospheres, 2nd ed., Freeman & Co., San Francisco 

Nakamura, T., Mazzali, P.A., Nomoto, K., & Iwamoto, K., 2001a, ApJ 550, 991 

Nakamura, T., Umeda, H., Iwamoto, K., Nomoto, K., Hashimoto, M., Hix, R.W., 
Thielemann, F.-K., 2001b, ApJ 555, in press 

Nomoto, K., Yamaoka, H., Pols, O.R., van den Heuvel, E.P.J., Iwamoto, K., Kumagai, S., 
& Shigeyama, T., 1994, Nature 371, 227 

Nomoto, K., et al. , 2001, in Supernovae and Gamma Ray Bursts, ed. M. Livio et al. 
(Cambridge: Cambridge University Press), in press ( |astro-ph/0003077| ) 

Patat, F. et al., 2001, ApJ 555, in press, ( |astro-ph/0103111| ) 

Richmond, M.W., et al. , 1996, AJ 111, 327 

Schlegel, E.M., & Kirshner, R.P, 1989, AJ 98, 577 

Sollerman, J., Kozma, C, Fransson, K., Leibundgut, B., Lundqvist, P., Ryde, F., & Woudt, 
P., 2000, ApJ 537, L127 



-20- 

Spyromilio, J., Stathakis, R.A., & Meurer, G.R., 1993, MNRAS 263, 530 
Spyromilio, J., 1994, MNRAS 266, 61 
Uomoto, A., 1986, ApJ 310, L35 



This manuscript was prepared with the AAS IATfjX macros v4.0. 



- 21 - 



Table 1. Parameters of the 'broad' synthetic spectra. 



date 


time from max 


SN epoch 


V 


M( 56 Ni) 


M(O) 


M(tot) 


D 7 t 


T e 


log n e 




days 


days 


km s _1 


M Q 


M Q 


M Q 




K 


g cm~ 3 


12 Aug 98 


+94 


108 


12500 


0.68 


3.0 


5.0 


0.31 


6286 


7.76 


12 Sept 98 


+125 


139 


12000 


0.62 


3.4 


5.2 


0.24 


5910 


7.47 


26 Nov 98 


+201 


214 


11000 


0.68 


5.0 


7.3 


0.19 


5092 


7.06 


12 Apr 99 


+337 


349 


10000 


0.65 


3.9 


5.4 


0.09 


4090 


6.40 


21 May 99 


+376 


388 


10000 


0.65 


3.2 


4.4 


0.07 


3874 


6.22 



t D 7 is the 7-ray deposition fraction 



Table 2. Parameters of the 'narrow' synthetic spectra. 



date 


time from max 


SN epoch 


V 


M( 56 Ni) 


M(O) 


M(tot) 


D 7 t 


T e 


log n e 




days 


days 


km s — 1 


M 


Mq 


Mq 




K 


g cm~ 3 


12 Sept 98 


+125 


139 


7500 


0.36 


2.1 


3.0 


0.32 


5971 


7.76 


26 Nov 98 


+201 


214 


6000 


0.33 


2.7 


3.8 


0.28 


5129 


7.43 


12 Apr 99 


+337 


349 


6000 


0.37 


2.0 


3.2 


0.12 


3987 


6.73 


21 May 99 


+376 


388 


5500 


0.34 


1.4 


2.2 


0.09 


3733 


6.61 



t D 7 is the 7-ray deposition fraction 
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Table 3. [O i] 6300 A fluxes and the O mass. 



SN epoch 


model 


F([0 I]) 


T 4 


M(0)(flux) 


M(O) (model) 


days 




erg cm -2 s _1 


10 4 K 


M Q 


M Q 


108 


broad 


6.0 10" 13 


0.63 


3.2 


3.0 


139 


narrow 


3.3 10~ 13 


0.60 


2.1 


2.1 


214 


narrow 


1.4 lO" 13 


0.51 


1.7 


2.7 


349 


narrow 


2.0 10~ 14 


0.40 


0.8 


2.0 


388 


narrow 


6.0 10~ 15 


0.37 


0.4 


1.4 
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i 1 1 1 1 1 1 1 1 1 1 1 r 




Rest Wavelength (A) 



Fig. 1.— The observed spectra of SN 1998bw from 12 August 1998 to 21 May 1999 compared 
to the 'broad lined' synthetic spectra (Table 1). All spectra have been shifted arbitrarily. 
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Fig. 2.— The observed spectra of SN 1998bw from 12 September 1998 to 21 May 1999 
compared to the 'narrow lined' synthetic spectra (Table 2). All spectra have been shifted 
arbitrarily. 



